makeblastdb -in ./R64.fa -input_type fasta -title SC_gDNA -dbtype nucl -out SC_gDNA
tblastn -query ./uniprotkb_taxonomy_id_4751_NOT_taxonomy_2024_03_13.fasta -db SC_gDNA -out ./blastx_reuslt.outfmt6 -evalue 1e-5 -outfmt 6 -max_target_seqs 1 -num_threads 11
./blast2gff.py -b blastx_reuslt.outfmt6 > blastx_reuslt.gff
sed  's/match_part/exon/g' blastx_reuslt.gff > blastx_reuslt_modified.gff
gffcompare -V -r R64.gff blastx_reuslt_modified.gff -o compare_result

